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ABSTRACT 



In this paper, we characterize the morphology of the disk-integrated phase 
functions of sateUites and rings around the giant planets of our Solar Sys- 
tem. We find that the s hape of the ph ase function is accurately represented 



5 by a logarithmic model (IBobrov 



197Ct in Surfaces and Interiors of Planets 



6 and Satellites, Academic, edited by A. Dollfus). For practica' 



param etrize the phase curves by a linear-exponential model (IKaasalainen et al. 



purposes, we also 



2001 



Journal of Quantitative Spectrosc opy and Radiative Transfe r, 70, 529 



19761 . Astronom- 



543) and a simple linear-by-parts model fiLumme and Irvind . 
ical Journal, 81, 865-893), which provides three morphological parameters : 
the amplitude A and the Half- Width at Half-Maximum (HWHM) of the oppo- 
sition surge, and the slope S of the linear part of the phase function at larger 
phase angles. 

Our analysis demonstrates that all of these morphological parameters are cor- 
related with the single scattering albedos of the surfaces. 
By taking more accurately into consideration the finite angular size of the Sun, 
we find that the Galilean, Saturnian, Uranian and Neptunian satellites have 
similar HWHMs (^0.5°), whereas they have a wide range of amplitudes A. 
The Moon has the largest HWHM (~2°). We interpret that as a consequence 
of the "solar size bias", via the finite size of the Sun which varies dramatically 
from the Earth to Neptune. By applying a new method that attempts to mor- 
phologically deconvolve the phase function to the solar angular size, we find 
that icy and young surfaces, with active resurfacing, have the smallest values 
of A and HWHM, whereas dark objects (and perhaps older surfaces) such as 
the Moon, Nereid and Saturn's C ring have the largest A and HWHM. 



26 Comparison between multiple objects also shows that Solar System objects 



2 



27 belonging to the same planet host have comparable opposition surges. This 

28 can be interpreted as a "planetary environmental effect" that acts to modify 

29 locally the regolith and the surface properties of objects which are in the same 

30 environment. 

31 keywords: Planetary rings; Satellites of Jupiter, Saturn, Uranus, Neptune, 

32 phase curves; opposition effect, coherent backscattering, shadowing, shadow- 

33 hiding, angular size of the solar radius 
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34 1 Introduction 



The opposition effect is a nonlinear increase of briglitness wlien tlie pliase an- 
gle a (the angle between the source of light and the observer as seen from 
the bod y) decrea s es to zero. This eff e ct wa s seen for the first time in Saturn's 



38 rings by 



Seeligerl (118841 ) and 



Miillerl (118851 ). Now, this photometric effect has 



39 been observed on ma ny surfaces in the Solar System : first on satellites of 



40 the gi ant planets, see ( 



41 oids, (iHarris et al 



Helfenstein et al. 



1989a 



42 Objects ( iBelskaya et al 



43 surfaces on Earth (IVerbiscer and Veverka 



1997) for a review; second on aster- 



Belskava and Shevchenko 



2008 



Rozenbush et al 



minerals in the laboratory ( 



Shkuratov et al 



2000l ) and Kuiper Belt 



20021): and finally o n various 



19901: IH 

L, 1999I 



Hapke et al 



Kaasalainen 



1996 )and for 



2OO3I ). The 



45 opposition effect on bodies in the Solar System has supplied interesting con 



46 straints about the regoli th and state of the surfaces (IHelfenstein et al. 



1997 



Mishchenko et al 



20061 ). Indeed, the opposition effect is now thought to be 
the combined effect of coherent backscatter (at very small phase angles), which 
is a constructive interference between grains with sizes near the wavelength 
of light, and shadow hiding (at l arger phase angle s ), wh ich involves shadows 



cast by the particles themselves (IHelfenstein et al. 



I997I ). 



By parametrizing the morphology of the phase functions for a ~ 0-20°, some 
numerical r nodels have derived physical properties of the medium in terms 



54 of regoliths (IMishchenko and Dlugachl. 



state of the macroscopic surface (IHapke 



1992a 



1986 



Shkuratov et al 



2002 



1999) and the 



Shkuratov et al. 



19991). 



However, such characterization of the phase function morphology is restricted 
by the angular resolution and the phase angle range of the observed phase 
function. Moreover, some effects (such as the finite size of the Sun and the 
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59 nature of the soil), which are not yet taken rigorously into account by the 

60 most recent models, can play important roles in a comparative study. 

61 For these reasons, it seemed important to test the behavior of the morphology 

62 of the phase function before using any physical model. 

63 The use of a simple morphological model is generally not adapted to derive 

64 the physical properties of the medium. But for the data set presented here, 

65 only the disk- integrated brightness I/F and the phase angle a are available, 

66 the corresponding angles of incidence (i) and angles of emission (e) are not 



67 given for these observations, so we cannot use 



tigations analytical models (IHapke 



1986 



s ophisticated for 



2002 



Shkuratov et al. 



urther inves- 



1999|) which 



need the brightness 1 jF and the three viewing geometry parameters a, [i and 
/io (/i and /io are the cosines of e and respectively). 

However, the theories developed for the coherent-backscattering and the shadow 
hiding effects deduce their properties b y parametrizing the opposition phase 



1999 



73 curve flMishchenko and Dlugach 



Hapke 



1986 



1992a 



Mishchenko 



1992 



Shkuratov et al. 



2OO2I ). Thus it is possible to connect the morphological pa- 



rameters A, HWHM and S with some physical characteristics of the medium 
derived from these models. 



77 The amplitude A of the opposition peak is general 

78 effects of the coherent-backscattering. According to 



y known to express th e 



Shkuratov et al. 



Nelson et al. 



fll999h 



(I2OOOI ). A is a function of grain size in such way that A decreases 
BO with increasing grain size (we refer to grains as the smallest scale of the sur- 

81 face compared to the wavelenght and virtual entities implied in the coherent 

82 backscatter effect, as microscopic roughness). This anti-correlation finds a nat- 

83 ural explanation in the fact that for a macroscopic surface, large irregularities 



84 with respect to the wavelength create less coherent effects than irregularities 
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85 with sizes comparable to t he wav elength. 



Mishchenko and Dlugachl (1l992bl ) and 



Mishchenkd (119921 ) emphasize that A 



B7 is linked to the intensity of the background A (defined as a morpholog ical pa- 



rameter of the linear-exponential function of 



Kaasalainen et al. 



89 is a decreasing function of increasing absorption (ILumme et al 



2001). which 



1990|); thus 



90 A must increase with increasing absorption or decr easing albedo vjo, w hich 



91 was c onfirmed by t he la boratory measurements of 



KaasalainenI (120031 ). In 



92 deed, iKaasalainenI (j2003l ) remarked that the opposition surge increases and 

93 sharpens when irregularities are small and that the opposition surge decreases 

94 with increasing sample albedo. 



The half width at half maximum HWHM, is also associated to the coherent- 
backscatter effect. It has been related to the grain size, inde x of refraction, and 



97 packing density of regolith, by p r evious numer ical studies (IMishchenko 



1992 



Mishchenko and Dlugach 



1992a 



Hapke 



2002). The varia tion of HWHM wit h 



19931 ) 



99 these three physical parameters is complex; see Fig. 9 of (IMishchenko 

100 HWHM reached its maximum for an effective grain size near A/2 and increases 

101 when the regolith grains' filling factor / increases. For high values of /, the 



103 However, several studies 



Helfenstein et al. 



(Il997h : 



2OO2I ) defined two HWHM parameters : for the 



Nelson et al. 


( 


2000) 


3 Hapke 


( 


2002 


) mo 



Hapke 



coherent- 
model of 



jackscatter HWHM (/ic), which is defined similarly to that in the 



Mishchenkd (119921 ). and the shadow hidi ng parameter hg. Applying 



107 this model to Saturn's rings. 



French et al. 



(120071 ) found that the coherent- 



backscatter peak is about ten times narrower than the shadow-hiding peak, 
but neither h^. nore hg equals to the morphological width of the peak HWHM. 
This reinforces the idea that a coupling of the two opposition effect mechanisms 
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Ill at small phase angles could be responsible for the observed surge width. 



112 Since the efficient regime of the shadow hiding is 



1985 



Helfenstein et al. 



1997 



Stankevich et al. 



0°-40°(lBuratti and Veverka 



19991) and that of the co - 



Helfenstein et al. 



19971) 



114 eherent bakscattering does not exceed several degrees ( 

115 the slope of the linear part S can be regarded as the only parameter that mir- 

116 rors the shadow hiding solely. This slope depends on the particle filling fac- 

117 tor D, which relates to the porosity of the regolith of a satellite and the ratio 



betw een the partic l e size and the phy s ical thickness of the ring 



119 ring (jlrvine 



1966 



Stankevich et al 



1999 



Kawata and Irvine 



or a p 



anetary 



19741 ). For a 



satellite, by "particles" we mean the macroscopic scales of the surface, which 
are implied in the shadow hiding effect . 



122 In the shadowing model of 



Irvind (119661 ) and 



Kawata and Irvind (119741 ) (which 



consists of the effects of shadows for a monolayer of particles), when the slope 
is shallow, the variation of a does not change the visibility of shadows and the 
particle filling factor must be high to make the proportion of shadows small 
for any observation geometry. By contrast, when S is steep, the particle filling 
factor is smaller and will contribute to a broad and large peak with a weak 
amplitude which will be regarded slope. 

In the shadow hiding model (i.e., multilayer shadowing), the larger the optical 
depth and the volume density D, the steeper the phase function is at large 



phase angles (10°-40°, 



Stankevich et al 



19991 ). However, at larger phase an- 



132 gles, the behavior of the absolute slope with albedo could change according 

133 to a more efficient regime of the shadow hiding (50°-90°, Stankevich 2008, 

134 private communication). 

135 For a compact medium such as a satellite's surface, the slope at very large an- 

136 gles (a>90°) is a consequence of topographic roughness, the so-called rough- 
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137 ness parameter 9 in the iHapkd (119841 . Il986l ) model. Then a steeper slope is 

138 due t o a surface tilt which varies from millimeter to centimeter scales (jHapkd . 



1984J ). However, the ro ughness can influence the ph ase curve at smaller phase 



angles as underlined by 



141 oratory measurements of 



Buratti and Veyerka fll985l ). also according to the lab- 



Kaasalainenl (120031 ) , the slope of the phase function 



142 (a<40°) increases with increasing roughness. From the theoretical assertions 

143 made above, the HWHM and the amplitude are governed by both coherent 

144 backscatter and shadow hiding effects, whereas the slope of the linear part of 

145 the phase curve is mirrors the unique expression of the shadow hiding effect. 

146 The goal of this paper is to understand the role played by the two known op- 

147 position effects (coherent backscatter and shadow hiding) on the morphology 

148 of the surge for different surface materials, which have different values of grain 

149 size, regolith grain filling factor, absorption factor (or inverse albedo), particle 

150 filling factor and vertical extension, by making some comparisons with the 

151 three morphological output parameters A, HWHM and S. 

152 This paper describes the results of a full morphological parametrization and 

153 comparison of phase functions of the main satellites and rings of the Solar 

154 System in order to compare the influence of parameters not yet implemented 

155 in actual models and simulations. Section 2 describes the data set that we 

156 used here and the specific reduction we added to these previously published 

157 data in order to compare them more easily. We also present the morphological 

158 models that have supplied the parameters and discuss their link with physical 

159 properties of the surfaces. In Section 3, we focus on the specific behaviors of 

160 the morphological parameters, as a function of the single-scattering albedo, 

161 the distance from the Sun and the distance from the center of the parent 

162 planet. Section 4 is dedicated to a discussion in which we physically interpret 
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163 the general behaviors obtained with a deconvolution method. Conclusions and 

164 future work that would be of interest are drawn in Sections 5. 

165 2 Data set description and reduction 

166 2.1 The opposition effect around a selection of rings and satellites of the 

167 giant planets 

168 We have applied a fitting procedure to a set of phase curves of satellites and 

169 rings obtained by previous ground-based and in situ optical observations (see 

170 Table 1 for references). The spectral resolution of the filters used for these 

171 observations are not rigorously mentionned by theirs authors, and because we 

172 mix for some objects phase curves of close wavelength, we give an approximate 

173 value of the wavelength of observation (the uncertainty of the approximated 

174 values is roughly 100 to 200 nanometers). 

175 For a comprehensive study of the morphology of the opposition phase curves, 

176 the solar phase curves of the Galilean satellites (lo, Europa, Ganymede and 

177 Callisto) and the jovian main ring were chosen, as well as the phase curves of 

178 the Saturnian rings (the classical A, B, and C rings and the tenuous E ring) 

179 and some Saturnian satellites (Enceladus, Rhea, lapetus and Phoebe); the 

180 rings and satellites of Uranus [We refer to the seven innermost satellites of 

181 Uranus - Bianca, Cressida, Desdemona, Juliet, Portia, Rosalind and Belinda 

182 - as the Portia group, to follow the designation of Karkoschka (2001). The 

183 phase function of the Portia group is then the averaged phase function for 

184 these seven satellites.], including the Portia group and three other Uranian 

185 satellites, Titania, Oberon and Miranda; and finally two Neptunian ring arcs 
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1B6 (Egalite and Fraternite) and two satellites of Neptune (Nereid and Triton). 

187 For all the satellites of this study, the phase function is representative of the 

188 leading side because they have, in general, better coverage at small phase 

189 angles (except for lapetus which have a trailing side brighter, we then use the 

190 lapetus' trailing side data). References for the phase curves that we use in this 

191 study are given in table 1. 

192 Insert Table 1 

193 This study should give an extensive comparison between rings around the 

194 giant planets (Jupiter, Saturn, Uranus and Neptune), as well as a comparison 

195 between rings and satellites for each giant planet of our Solar System. For 

196 practical purposes, the well-known phase curve of the Moon is added as a 

197 reference. 



198 2.2 Data set reduction 



199 In order to properly compare the morphological parameters of the objects 

200 whose phase curves are given as magnitudes, we have converted the magni- 

201 tude M to the disk-integrated brighness I/F hj using : 



I/F = 10 



-OAM 



202 ( jPomingue et al. 



19951 ). This modification allows us to directly compare the 



203 slope of the linear part of all the curves in the same unit. 
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204 2.3 Data set fits: the morphological models 



The purpose of the present paper is to provide an accurate description of the 
morphological behavior of the observed phase curves. This is the very first step 
prior to any attempt to perform either analytical or numerical modeling. As 
a consequence, special care has been given here to parametrizing the observa- 
tions efficiently and conveniently. In addition, morphological parametrization 
is necessary to efficiently compare numerous phase curves and derive statisti- 
cal behavior, as will be done in Sections. 

Several morphological models have been used in the past to quanti tatively 



213 descri be the shape of the phase fun ctions : the logarithmic m odel of 



Bobrov 



214 ( 1l970l ). the linear- by- 



215 exponential model of 



parts model of 



Kaasalainen et al. 



Lurnme and Irvind (119761 ) and the linear 



(120011 ). The specific properties of these 



three models make them adapted for different and complementary purposes. 
The logarithmic model is an appropriate and simple representation of the data, 
the linear-by-parts model is convenient to describe the shape in an intuitive 
way, and finally the linear-exponential model is commonly used for the phase 



220 curve s of Solar System bodies (see the comparative study of 



Kaasalainen et al. 



20011 ). 



222 Insert Fig. 1 



223 2.3.1 The linear-by-parts model 

224 For an intuitive description of the main features of the phase curves, the linear- 

225 by-parts model is the most convenient one. It is constituted of two linear 

226 functions fitting both the surge at small phase angles (a < ai) and the linear 
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227 regime at larger phase angles (a > 02), where, generally, ai 7^ a2- Besides ai 

228 and a2, this function depends on 4 parameters, Ao, Bo, Ai, and Bi, such that: 



I/F{a < ai) 
I/F{a > 02) 



-Ao ■ a + Bq 
-Ai - a + Bi 



(2) 
(3) 



Lumme and Irvind (119761 ) and 



Esposito et al. 



(11979h use ai=0.27° and a2=1.5°. 



230 By testing several values of ai, it appears that for our data set, values of 

231 ai = 0.5° and 0:2=2° provide the best results, so these values are now adopted 

232 in the rest of the paper except for the Moon and tenuous rings for which we 

233 take ai = 1°. 

234 In terms of the four parameters Aq, Bq, Ai, and Bi, the shape of the curve 

235 is characterized by introducing three morphological parameters : A, HWHM 

236 and S designating the amplitude of the surge, the half-width at half-maximum 

237 of the surge, and the absolute slope at "large" phase angles (i.e., a few degrees 

238 up to tens of degrees) , respectively. The parameters are defined by : 



A=^ HWHM= ^!^rV. and S = A, (4) 
Bi 2{Ao -A,) ^ ^ 

239 Even if all the part of the opposition curve cannot be fitted by two linear 

240 functions, this model offer a convenient description of the main trends of the 

241 phase curve. 



242 2.3.2 The linear- exponential model 

243 The linear-exponential model describes the shape of the phase function as a 

244 combination of an exponential peak and a linear part. Its main interest is that 

245 it has been used in previous work for the study of the backscattering part of the 
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246 phase curves of the Solar System's icy satellites and rings (IKaasalainen et al. 



2001 



Poulet et al. 



2002h. 



248 However, as noted by 



French et al. 



(120071 ). we find that this model does not 



fit the phase curves well: in particular A, HWHM and S are under- or over- 
estimated. In addition, the converging solutions found by a downhill simplex 
technique have large error bars, which means that a large set of solutions is 
possible and thus produce some difficulties for the comparison with the other 
objects. 

For completeness, we give the four parameters of this model : the intensity 
of the peak /p, the intensity of the background h, the slope of the linear 
part Is and the angular width of the peak w such that the phase function is 
represented by : 



I/F = h + Is-a + Ip- e-^ (5) 

258 As a — > 0, exp{—a/2w) ^ 1 — a/2w + C(a^), so that the slope approaches 

259 Is + Ip/{2w). The degeneracy of these parameters may explain some of the 

260 difficulty in obtaining good fits described above. For consistency with previous 

261 work, we can express the amplitude and HWHM of the opposition surge in 

262 this model as: 



^^Ip+Jb HWHM = 2-ln2w; and S = -h (6) 



263 We report in Table 2, the morphological parameters of the linear-by-parts 

264 model and that of the linear-exponential model. 

265 Insert Table 2 
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266 2.3.3 The logarithmic model 



267 As no ted by 



Bobrovl fll970h 



Lumme and Irvind (119761 ) and 



Esposito et al 



268 ( 1l979l ). we remark that a logarithmic model describes the phase curves very 

269 well. It depends on two parameters (oq and ai). This model has the following 

270 form : 



I/F = ao + ai ■\n{a) (7) 

271 In general, this model is the best morphological fit to the data. However, 

272 oo and ai are not easily expressed in terms of A, HWHM and S, since the 

273 model's dependence on a is scale-free. Thus we report the values of these two 

274 parameters in table 3 to allow an easier reproduction of the observational data. 



275 2.3.4 ^ method that takes into account the angular size of the Sun 

276 For all the phase curves presented here, a comparison of their surges could be 

277 compromised because they have different values of their observed minimum 

278 phase angle values. For example, data for the Galilean satellites never reach 

279 0.1°, whereas data for Saturn's rings almost reach 0.01°. 

280 Although the behavior within the angular radius of the Sun represents a small 

281 part of the phase function, these smallest phase angles are crucial to constrain 

282 the fit, especially for the linear-exponential model. Indeed, when a ^ 0, the 

283 linear-exponential function tends toward + Ip/{2w)] as a consequence, this 

284 function flattens at very small phase angles. 

285 However, in some cases this flattening does not correspond to the expected 

286 flattening due to the angular size of the Sun because the linear-exponential 
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2B7 flattening fits itself arbitrarily with the phase angle coverage. The less points 

288 there are at small phase angles, the sooner will occur the flattening of the 

289 phase funct i on. 



DeauetaL 



fl2008l ) showed for Saturn's rings that the behavior of the surge was 

291 accurately represented by a logarithmic model between 15° and 0.029°, where 

292 0.029° corresponds to the angular size of the Sun at the time of the Cassini 

293 observations. Below 0.029°, the resulting phase function flattens, whereas the 

294 logarithmic function continues increasing. The use of this fact observed specif- 

295 ically for Saturn's rings and its generalization to the Solar System objects of 

296 this study allows us to create extrapolated data points. Indeed, it is more con- 

297 venient to use extrapolated data points than convolve the linear-exponential 

298 function or the logarithmic function with the solar limb darkening. First be- 

299 cause, for incomplete phase functions, the flattening of linear-exponential func- 

300 tion is almost uncontrollable and second because for the logarithmic model, 

301 even if a convolution is possible, linking the morphological parameters A, 

302 HWHM and S to the outputs Qq and ai is not trivial. 

303 The method to create extrapolated data points consists of flrst fltting the 

304 logarithmic model to the data and then taking the value of the logarithmic 

305 function at the phase angle which corresponds to the solar angular size (a©, 

306 see Appendix). We then give to six points the same y- value : l/F{a = a©) 

307 and x-values ranging from 0.001° to of phase angle. These extrapolated 

308 data points are represented in Figure 2 (the full method is detailed in the Ap- 

309 pendix). The extrapolated data and the original data are then fltted by the 

310 linear-exponential model in the last step. 

311 Insert Fig. 2 

312 For the Moon, Ariel and Oberon, for which the phase curve has a few points 
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313 below the solar angular radius, we can see that the extrapolated points match 

314 the observational points quite well. This proves that the solar angular size 

315 effect is a ffattening of the phase function below a©. In the case of the HST 

316 data for Saturn's rings, for which we have also a few points below the solar 

Dserved 



French et al. 



fl2007f ) 



317 angular radius, the extrapolated points are a bit smalle r than the o 

318 points. This is may be due to the fact that the data of 

319 are already deconvolved by another method. 

320 We also performed a convolution of the linear-exponential function to a limb 

321 darkening function (see Appendix), but this refinement did not significantly 

322 change the values of A, HWHM and S, because the linear-exponential model 

323 already flattens as a ^ 0. Thus, by adding extrapolated data below a©, we are 

324 sure that the resulting fitting function will have a constant behavior below a© 

325 and that the resulting fitting function will take into account the angular size 

326 of the Sun. However, we assume that all bodies have a logarithmic increase up 

327 to the solar angular radius, which is only confirmed for Saturn's rings. Output 

328 parameters of our best fit for the "extrapolated linear-exponential" function 

329 are given in tables. 

330 Insert Table 3 



331 2.3.5 A method of solar size deconvolution 

332 Although the behavior at phase angles smaller than the solar angular radius 

333 represents a small part of the surge, a comparison of the surge of Solar Sys- 

334 tem rings and satellites could be compromised because they have different 

335 values for the mean solar angular radius (anim=0.051, 0.028, 0.014, 0.009° re- 

336 spectively for Jupiter, Saturn, Uranus and Neptune at their mean distances 
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337 from the Sun). Indeed, according to the results presented here, the amplitude 

338 and HWHM seems linked to the finite angular size of the Sun. However, our 

339 morphological study doesn't clearly show that the Sun's angular size effect is 

340 preponderant for the amplitude A, because even considering more accurately 

341 the "solar environmental effect," the surges of Neptune's satellites have smaller 

342 amplitudes than those of Uranus (figure 8b). This contradicts the theoretical 

343 assumption that the solar angular size would give the largest amplitude to the 

344 most distant objects, for which the Sun has the smallest angular size. Because 

345 we previously noted that the effect of the "solar environmental effect" was to 

346 fiatten the phase fun ction when the p hase angle is less than or equal to the 



347 solar angular radius (jPeau et al. 



20081 ). a naive deconvolution method would 



348 be to allow the phase curve to rise below a©. This is also suggested by a pre- 

349 vious deconvolution of HST data on S aturn 's rings, for which the brightness 



350 still increases below ap, ( French et al. 



20071 ). However, the linear-exponential 



351 function is not appropriate for this purpose because it intrinsically flattens as 

352 a — > 0. Thus the only morphological function that allows an increase, even at 

353 very small phase angles, is the logarithmic function. In particular, this function 

354 allows the same increase of the brightness above and below the solar angular 

355 radius (without break in the brightness), then using this function simulates a 

356 point source of light. However, we assume that the smallest phase angles are 

357 about a=0.001°, in this way, if a physical flattening should be performed by 

358 the coherent backscattering or the shadow hiding effect, it will be possible at 

359 these phase angles. 



360 The logarithmic functi on flts the phase 



361 phase angles (0.1-15°, 



Deau et al. 



2008 



unction quite well at small and large 
); however, none of the morphologi- 



362 cal parameters A, HWHM and S are well-deflned in this model. Because the 

363 logarithmic model is a good representation of the data and perform a kind 
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364 of deconvolution at phase angles less than Oq, we fit this function by the 

365 linear-by-parts function. As shown in figure 3, the fitting results of this crude 

366 deconvolution method is reasonably acceptable when the phase function is 

367 plotted on a linear scale of phase angle. 

368 Insert Fig. 3 

369 However, on a logarithmic scale of a, the linear-by-parts fit is obviously not 

370 acceptable because intrinsically, a linear function cannot fit a logarithmic in- 

371 crease. As a consequence, we have slightly changed the linear-by-parts pa- 

372 rameters in order to take into account the inappropriate flattening of this 

373 function, compared to the logarithmic function. Because the y-intercept Bq of 

374 the linear-by-parts model is less than the values of the logarithmic function 

375 when a<0.01°, we replace Bq by the value of the logarithmic function when 
a=0.001° : 



5^ = ao + ai-ln(10~^), (8) 

377 where oi < 0, in such a way that the amplitude and the angular width are 

378 now given by : 



A=^ and HWHM = 4^V-^ (9) 
Bi 2{\ai\-Ai) ^ ' 

379 where |ai| is the absolute slope of the logarithmic function (see section 2.3.3). 

380 Replacing |ai | by Aq in the formula for HWHM was motivated by the fact that 

381 the original values of the linear-by-parts parameters were systematically the 

382 same (HWHM~0.22° for ai=0.3°). This is due to the fact that the logarithmic 
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function is a fractal function, so it is not possible to obtain a Half- Width at 
Half Maximum. Values of the linear-by-parts parameters A, HWHM and S 
are given in Table 5. Our best result is probably th at for the B ring, for which 



we previousl y found different va. 
data and the 



French et al 



ues of A from the 



Franklin and Cookl (119651 ) 



(120071 ) data (A=1.30 and A=1.38 respectively, see 



table 2 with the convolved models). The discrepanc y was still present w i th the 



3B9 extrapolated lin ear-exponenti a 



and A=1.32 for 



French et al. 



model (A=1.35 for 



Franklin and Cookl (119651 ) 



(120071 )). due to the fact that the solar angular 
size was different at the two observation times (see table Appendix6). Now, 
with the unconvolved model, the discrepancy of the two values is somewhat 



393 reduced compared to va. 



Franklin and Cookl (119651 ) and A=1.77 for 



ues from the fit to the original dat a : A=1.82 for 



French et al 



(120071 ). table 5, which 



395 implies that the angular size effect is now absent. 



396 Insert Table 5 



397 3 Results 



398 Our procedure is to interpret more carefully the morphological results of a 

399 large dataset. We start first by studying the behaviors of the morphological 

400 parameters with the single scattering albedo with the raw data (Section 3.1 

401 and 3.2) and the improved data that take into account of the solar size of the 

402 Sun (Section 3.3). In a last step, we freed from the solar size biais by trying to 

403 look the opposition effect in the outer Solar System with the same solar size, 

404 assumed to be a point (Section 3.4). 
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405 3. 1 Behaviors of the morphological parameters 



406 In this section we compare the morphological parameters as a function of the 

407 single scattering albedo Wq. 



Since the single scattering albedo represents the ratio o: 



409 ciency to total light extinction over the phase angle range (j Chandr asekhar 



scattering effi- 



410 Il960l ). its value must be computed with the largest coverage of phase angle as 

411 possible (0 to 180 degrees). This is why we did not compute the single scatter- 

412 ing albedo with the phase curves presented in this paper but we use previously 

413 published values of single scattering albedo computed from phase curves with 

414 a larger phase angle coverage than ours and a wavelength close to ours. Thus, 

415 references for phase curves (table 1) and references for wq (table 4) are not 

416 always the same. 

417 Insert Table 4 

418 We did not found single scattering albedo values for the jovian main ring and 

419 the Saturn's E ring, so these two objects will be excluded of the study of the 

420 mophological parameters with the single scatterig albedo. 



421 3.1.1 Variation of the angular width of the surge with albedo 

422 First, we discuss the variation of HWHM=/(wo) derived from the linear-by- 

423 parts model (Figure 4a) and HWHM=/(roo) derived from the extrapolated 

424 linear-exponential model (Figure 4b). Interestingly, the variation differs ac- 

425 cording to the morphological model : the first case leads to a decrease of 

426 HWHM when wq increases, while the latter case leads to an increase of HWHM 

427 when wq increases. 
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428 Insert Fig. 4 



429 The different results for HWHM=/(ti7o) between the linear-by-parts results 

430 (well fitted by HWHM~ 0.9 x 0.3"^°, Figure 4a) and the extrapolated linear- 

431 exponential results (represented by HWHM~ 0.15 + O.ldwo, Figure 4b) is 

432 mainly due to the points which correspond to the Moon, Callisto and Nereid. 

433 Indeed, in general, values from the extrapolated linear-exponential model sig- 

434 nificantly decrease for the outer Solar System objects, whereas the value for 

435 the Moon increases by almost 1°. This is due to the fact that when we take 

436 into account the Sun's angular size, this effect lower the values of HWHM for 

437 the incomplete phase functions. 

438 However, Figure 4 shows a large dispersion of HWHM with albedo. 



439 3.1.2 Variation of the amplitude of the surge with albedo 



440 Figure 5 shows a weak dependence of the amplitude of the s u rge on the albedo 



441 for th e satellites, already noted by 



Helfenstein et al. 



(120021). 



(11997h 



Rozenbush et al. 



443 Insert Fig. 5 



444 In both cases (linear-by-parts model. Figure 5a and extrapolated linear-exponential 

445 model. Figure 5b) we note a dependence of A with O7o, which follow a function 

446 leading to a decrease of A when vuq increases (A~ 1.65 x 0.72™" in Figure 5a 

447 and A~ 1.75 x 0.72^° in Figure 5b). The consistent trends in both cases im- 

448 ply that the finite size of the Sun was correctly derived by the linear-by-parts 

449 model. 

450 The decrease of A with increasing wq could be understood by a relation be- 
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tween the amplitude and the single scattering albedo via the intensity of the 



453 verse 



) , which is in- 



Lumme et al. 



background phase function lb (with the linear-exponential mode! 
y proportional to the albedo. Thus, the predicted trend of 

454 (jl990l ) is confirmed by our present results. 

455 In addition, Figure 5b labels satellites by color to indicate their parent planet. 

456 This figure indicates that distant objects (such as the Uranian satellites) have 

457 a significantly larger amplitude than less distant objects (such as the Galilean 

458 or Saturnian satellites) while the Neptunian satellites have values in the av- 

459 erage. Then it must be considered that the finite size of the Sun has a role in 



the amplitude's value fIShkuratov 



199ll ). As a consequence, even if the trend 



of A=/(ti7o) is well explained by theoretical considerations, one can remark 
that the large dispersion in this correlation could be due to other effects (such 
as the finite size of the Sun) that weakens the albedo dependence of A. Thus, 
as for HWHM, we cannot physically interpret the variation of HWHM and A 
as long as they are convolved with the effect of the solar angular size. 



466 3.1.3 Variation of the slope of the linear part with albedo 

467 The last morphological parameter is the slope S, which we represent as a func- 

468 tion of the single scattering albedo for the rings (Figure 6a) and satellites 

469 (Figure 6b) of the Solar System. 

470 Insert Fig. 6 

471 In this figure, rings and satellites have different values of slope as function 

472 of their albedo, and a slight increase for S with increasing vuq is noticed. For 

473 the rings (Figure 6a), it seems that a good correlation appears between S and 

474 the albedo, which may be roughly fitted by a function like S~0. 001+0. 02-'a7o- 
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475 A similar fit works well for the satellites (the Moon, Saturnian and Uranian 

476 satellites are not far from the dashed line in figure 6b). This fit to the points 

477 could be S~0. 001+0. OI-cctq (Figure 6b). However, three objects fall far from 

478 this curve : Europa, Ganymede and lo. This correlation suggests that multiple 

479 scattering may be a strong element at play in the regim e of self-shadowing 



(beyon d ~1° of phase angle), in qualitative agreement with 



Kawata and Irvine 



fll974h . 



4B2 3.2 Cross comparisons between the morphological parameters of the surge 

4B3 We see in Figure 4 that the angular width of the surge can change significantly 

4B4 by taking into account the solar angular radius. However, it is not the case for 

485 the amplitude of the surge (Figure 5). Figure 7 shows the behavior of a cross 

486 comparison between the morphological parameters A and HWHM obtained 

487 with the linear-exponential model convolved (Figure 7b) or not (Figure 7a) 

488 with the limb darkening function. 

489 Insert Fig. 7 

490 In the first graph (Figure 7a), it first seems that two different groups may be 

491 qualitatively distinguished. 

492 On the one hand, there is a group of objects with similar values of the HWHM, 

493 in the range 0.1° to 0.4°, but with significantly different values of the ampli- 

494 tude, from 1.4 to 1.8. It is interesting to note that these bodies, which include 

495 the Saturnian rings and Uranian satellites, are not bodies in the outermost 

496 part of the Solar System (such as the Neptunian satellites). Within this group, 

497 we also note that similar objects are gathered in the (A, HWHM) space: the 
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498 Uranian satellites have, on average, the largest values of the amplitude, ~1.7. 

499 Saturn's rings have an amplitude between 1.3 and 1.6, closer to the Uranian 

500 satellites. We also note that whereas all satellites have quite a constant HWHM 

501 (between 0.2° and 0.4°), Saturn's rings have systematically lower values, be- 

502 tween 0.08° and 0.09°, which may be suggestive of a different state of their 

503 surface. 

504 The second group includes Saturn's satellites, along with lo, Europa and Tri- 

505 ton, which have the lowest values of amplitude. A very striking feature is the 

506 peculiar behavior of bodies such as Callisto and the Moon : they have similar 

507 amplitudes (about 1.5) and also similar HWHMs (about 2°). In the second 

508 graph (Figure 7b), it first seems that the bodies belonging to the same pri- 

509 mary planet have similar values of A and HWHM. For the Galilean satellites, 

510 we found the largest HWHM for the outer Solar System satellites (between 

511 0.2° and 0.5°) and amplitude between 1.1 and 1.5. The Uranian satellites still 

512 have the largest amplitudes (A ranges between 1.6 and 1.9), but the values 

513 of HWHM are similar to that of those of Saturn's rings and satellites. The 

514 Neptunian satellites have the sharpest opposition peaks (HWHM<0.1°) but 

515 moderate amplitudes (between 1.2 and 1.4), similar to the range of the Sat- 

516 urnian satellites. 

517 Does this imply some deep structural difference of the surface regolith of bod- 

518 ies, or is it due to the Sun's angular size effect? For the moment we note that 

519 the opposition effect is poorly understood, especially at phase angles smaller 

520 than 1° in the coherent backscattering regime. 

521 Whereas physical implications are still hard to draw from these graphs, it is 

522 interesting to note that the solar angular size refinement that we use naturally 

523 clusters different kind of surfaces in different locations of the (A, HWHM) 

524 space, and that "endogenically linked objects" are quite well gathered in small 
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525 portions of this space. This could suggest that common environmental pro- 

526 cesses (meteoroid bombardment, surface collisions, space weathering, etc.) may 

527 homogenize different surface states by processing mechanisms that may deter- 

528 mine the microstructure of the surface, and then, in turn, the behavior of 

529 the opposition surge at very low phase angles, as it may be linked wi t h the 



530 spatial organization of m i crometer-scale surface regolith (IMishchenko 



1992 



Mishchenko and Dlugach 



1992a : 



Shkuratov et al. 



19991). 



532 3.3 Additionnal effect 



533 With the cross comparison of the morphological parameters of the surge, the 

534 angular size of the Sun and the fact that objects seem "endogenically linked", 

535 two supplementary effects (observationnal and physical) can significantly mod- 

536 ify the values of A and HWHM : the "solar size bias" and the "planetary 

537 environmental effect". 



538 3.3.1 The "solar size bias" 

539 We tested the influence of the solar angular size by representing in Figure 8 

540 the morphological parameters of the surge A and HWHM as a function of the 

541 distance from the Sun d (in Astronomical Units). 

542 Insert Fig. 8 

543 We represent HWHM (Figure 8a) and A (Figure 8b) from the linear-by-parts 

544 model with filled symbols and that of the extrapolated linear-exponential 

545 model with empty symbols. 

546 We remark that the linear-by-parts angular width follows the power-law func- 
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547 tion HWHM~0.33+l.l(i-i-5 (the solid line in Figure 8a). The fit is quite good 

548 from the Moon to Uranus, but is far from the values of Neptune's satellites 

549 (especially that of Nereid). It is easier to see with this representation that 

550 the HWHM of Nereid is larger than that expected by the power-law function. 

551 The extrapolated linear-exponential HWHM corrects this because the Nep- 

552 tunian satellites now have smaller values of HWHM that are better fitted by 

553 the power-law function. Indeed, the extrapolated linear-exponential HWHM 

554 follows a similar function (HWHM~0.12+2.3(i~^'^), but values at the extreme 

555 parts of the Solar System (innermost with the Earth's satellite and outermost 

556 with Neptune's satellites) are significantly different : for the Moon, the extrap- 

557 olated linear-exponential HWHM is larger than its linear-by-parts counterpart 

558 and for the Nereid, the extrapolated linear-exponential HWHM is smaller than 

559 the linear-by-parts HWHM. 

560 However, such a strong trend is not observed in the case of the amplitude 

561 of the surge. As shown in Figure 8b, a fit to the linear-by-parts amplitudes 

562 is good for the Galilean, Saturnian and Uranian satellites (which we fit by a 

563 linear function A~l.l+0.08(i) but not at all for the Moon and the Neptunian 

564 satellites. The predicted behavior (dashed line in Figure 8b) shows that the 

565 value for the Moon is overestimated and that the values of the Neptunian 

566 satellites are strongly underestimated. The use of values of the extrapolated 

567 linear-exponential A did not improve the fit. Indeed, the extrapolated linear- 

568 exponential amplitude is larger than the linear-by-parts amplitude for the 

569 Moon, whereas the extrapolated linear-exponential values of the Neptunian 

570 satellites are smaller than their linear-by-parts counterparts. The exact op- 

571 position trends were expected to obtain a good linear fit from the Moon to 

572 Neptune. Perhaps the solar size effect is not important at Neptune's distance 

573 and the values of A and HWHM are physical, in the sense that they depend 
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574 only on the opposition effect mechanisms (coherent backscattering and shadow 

575 hiding) . 

576 These results seems to suggest that A is less affected by the "solar size bias" 

577 than HWHM, which is entirely controlled by this effect (which seems trivial 

578 because this bias is an angular effect). It is possible that the values of A result 

579 from a coupling of the physical opposition effects (coherent backscatter and 
5B0 shadow hiding) with the environmental opposition effects (solar and plane- 

581 tary). As a consequence, the deconvolution of the phase function (at least for 

582 the "solar size bias") should allow the physical opposition effects to express 

583 fully themselves in the values of A and HWHM. 



584 3.3.2 The "planetary environmental effect" 



585 Previous studies by 



Bauer et al. 



mm and 



Verbiscer et al. 



(120071 ) have con- 



586 firmed, at the scale of the Saturnian system, a kind of "endogenic" or "ecosys- 

587 temic" classification of the opposition surge. Indeed, these works demonstrated 

588 that the opposition surge paramaters of the outermost and innermost Satur- 

589 nian satellites, respectively, can be a function of the distance from Saturn. 

590 For the planetary environments of Jupiter, Uranus and Neptune, there is no 

591 significant variations with distance from the parent planet. The first reason is 

592 maybe statistical because there are not enough data to make (for these sys- 

593 tems, we have less than four objects). The second is that the dust environment 

594 can be influenced by other effects (the magnetospheric activity, the satellite's 

595 activity, the proximity to the Kuiper Belt and transneptunian objects), in such 

596 way that the distance from the planet host could be irrelevant for some of the 

597 planetary environments. 

598 For the Saturnian system, for which local interactions between satellites and 
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599 rinKs exists^ we ob serv e similar trends with our results than the trends of 



Bauer et al. 



(120061) and 



Verbiscer et al. 



(120071) (see Figure 9). 



Insert Fig. 9 

Variations of the morphological parameters on large scales of distance (for the 
Saturn system) show trends that suggest a common ground for environmental 
processes. These processes may imply different surfaces, but will be handled 
by the opposition effect in the same way by the mechanisms that determine 
the microstructure of the surface. According to theoretical models of coherent 
backscattering, the amplitude is related to the grain size and HWHM depends 
on the composition, distribution of grain size and the regolith filling factor. 



Thus the behavior o 



the opposition surge is connected to the spatial organiza - 



610 tion of the regolith (IMishchenko and Dlugach 



1992a 



Shkuratov et al. 



19991). 



Therefore, the study of the morphology of the opposition peak can highlight 
dynamical interactions between the rings, satellites and the surrounding en- 
vironment through the photometry. These ring/satellite interactions noticed 
here go beyond the general dynamical interactions between rings and satel- 
lites (such as resonances, for example). Here these interactions involve common 
erosion histories on the surfaces of the rings and satellites. S i milar v alues of 



617 HWHM according to the theory of 



Mishchenko and Dlugachl (jl992al ) can be 



explained by similar values of refractive index (with various values of grain 
size and filling factor), or by different values of refractive indices, but similar 
values of grain sizes and filling factor of the regolith. 

There are two known mechanisms that can act together in order to explain the 
similarities in the values of HWHM for objects that are endogenically linked. 
The impacts of debris in planetary environments can change the chem- 
ical composition of the rings and satellites : new elements can be directly 
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625 added to the system ; the more volatile elements can be preferentially re- 



626 moved and the more fragi 



work of 



e comp ounds can be preferentially processed. The 



Cuzzi and Estrada] (119981 ). in particular, details changes in the chem- 



ical composition of Saturn's rings by meteoroid bombardment and ballistic 
transport. 

(2) The second mechanism that is likely to act concerns every kind of coUi- 
sional mechanism capable of modifying, at microscopic scales, the surface of 
the satellite's regolith (met eoroid bombardment , external collisions, d i sinte 



633 grati on in space, etc.); see (ILissauer et al. 



1988 



Colwell and Esposito 



1992 



19931 ). In the case of ring particles, models of erosion by ballistic transport 



635 have been developed an d predict the destruction of micrometer-sized grains in 



dense rings (r> 1), see (lip 



19831). 



637 3.4 Study of the unconvolved opposition parameters 

638 With the unconvolved morphological parameters obtained with the method 

639 in Section 2.3.5, we are now sure that the morphological surge parameters are 

640 independent of the distance from the Sun, and thus independent of the "solar 

641 size bias". Indeed, in Figure 10, we represent A and HWHM derived from the 

642 linear-by-parts model which fits the logarithmic model as a function of the 

643 distance from the Sun and there is no relation with the distance, unlike in 

644 Figure 8. 

645 Insert Fig. 10 

646 Moreover, we can see that three groups can be distinguished : 

647(1) A group with the Moon, Callisto, the C ring, Phoebe, the classical Uranian 
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satellites (Ariel, Titania and Oberon) and Nereid. They have large angular 
widths and amplitudes (HWHM>0.5°, Figure 10a and A>2.0, figure 10b). 
We can see that these objects are dark, with low and moder ate albedos (ta- 



ble4). These objects are also known to be heavily cratered (INeukum et al. 



2001 



Zahnle et al. 



20031 ): thus, they don't have intrinsic resurfacing mech- 
anisms. For the specific case of the resurfacing of the Saturn's C ring, it 
is known th at the collisional activity of a ring is controlled by the opti- 



655 cal depth r (ICuzzi and Estrada 



19981 ). The number of collisions per orbit 



656 per particle is proport i onal 



Wisdom and Tremaine 



to r (in the regime of low optical depth, see 



19881 ). and the random velocity in a ring of thick- 

658 ness H is about H x Q (with Q standing for the local orbital frequency). 

659 Since H is a decreasing function of r, impact velocities are high in regions 

660 of low optical depth. As a result, particles in low optical depth regions (such 

661 as the C ring) may suffer of resurfacing characterized by rare, but somewhat 

662 higher-speed, collisions. 

663(2) A group with lo, lapetus, Rhea and the bright Saturnian rings character- 

664 ized by smaller amplitude and angular width : A<1.7 and HWHM~0.4° 

665 (figures 10a,b). 

666 ( 3) A group with Ganymede, Europa, Enceladus and Triton with the smallest 

667 amplitude and angular width (1.3<A<1.6, Figure 10a and 0.1<HWHM<0.3°, 

668 Figure 10b). Interestingly for the amplitude, we can see that this group con- 

669 tains only the brightest surfaces of the Solar System, with a single scattering 



albedo close to zuq_^0.9 



see 



howey er, Bond albedo of Ganymede is quite lower. 



Squyres and Veverkal . 



198ll ). These objects are also known to have active 



672 resurfacing. Indeed, this was confirmed for Europa, which h as a very young 



673 surface and perhaps recent gey ser-like or volcanic activity (ISullivan et al. 



1998 



Pappalardo et al 



1998bl ). and Ganymede, on which the grooved ter- 
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675 rains co uld have formed th r ouKh tectonism, probably co mbined with icy vol- 



676 canism (jPappalardo et al. 



1998 



McCord et al. 



200ll ). Present-day resur- 



677 facing is also taking place on Enceladus, whose geysers produc e the E ring 



For CO et al 



19951 ) ■ 



20061 ). and for Triton, which also has geysers fl Croft et al. 



680 This classification seems to suggest that the darkest and oldest surfaces have 

681 the largest amplitudes for the surge and that the brightest and youngest sur- 

682 faces have the smallest amplitudes. However, one might be surprised that the 

683 third group does not include lo, which has an intense resurfacing via tidally 

684 induced volcanism. Also, the Portia group does not belong to only one group 

685 in figure 10 : it belongs to the group 3 for HWHM and to the group 2 for 

686 the amplitude A. For these two isolated cases, it is possible that the average 

687 of photometric phase curves from different satellites is responsible of the fact 

688 that lo and the Portia group are difficult to classify. 



689 4 Discussion 



690 4-1 Implications of the surge parameters of the unconvolved data 

691 By removing the "solar size bias", we can try to physically interpret the am- 

692 plitude variations with the single scattering alebdo with the mechanisms pro- 

693 posed to explain the opposition effect. Our study shows a link between the 

694 single scattering albedo and the unconvolved morphological parameters. A 

695 linear fit to the unconvolved amplitude is A=2.2 — 0.5ti7o (with a correlation 

696 coefficient of -61%) and the linear fit to the unconvolved angular width is 

697 HWHM=0.52 — 0.19tao (with a correlation coefficient of -38%). By excluding 
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the Portia group, we find a better correlation coefficient for HWHM : -66%. 
These correlations are stronger than that previously found with the convolved 
data. This shows that the "solar size bias" acts to scatter the morphologi- 
cal parameters. As a consequence, the fact that old and dark surfaces with a 
low resurfacing activity have high unconvolved amplitude whereas the bright 
and young surfaces with an intense resurfacing activity have low unconvolved 
amplitude is linked to the single scattering albedo variations of A. Indeed 



the single scattering a 
not only: according to 



bedo is a measu r e of t he brightness of a surface. But 



Shkuratov et al. 



(Il999l ). the amplitude of the coherent 



backscattering opposition surge is a decreasing function of increasing regolith 
grain s ize. If the morphological amplitude is due to the coherent ba ckscattering 



709 effect f lMishchenko and Dlugach 



1992b 



Mishchenko et al 



20061 ). the depen- 



dence of A = fiwo) could be understood as a positive correlation between 
the grain size and the single scattering albedo. However, it is possible that the 
morphological amplitude is not only that of the coherent backscattering effect 



713 but is domina ted by 



714 underlined by 



3 oth ef fects : coherent backscatter and shadow hiding, as 



Hapkd (I2OO2I ). However, here it is not possible to separate the 



two effects and say which effect is dominant because to separate the coher- 
ent backscatter and sh adow hiding mechanisms, the polarization is required 



( iMuinonen et al. 



20071). 



718 4-'^ Implications of the slope of linear part 

719 The strong correlation of the slope S (in I/F.deg~^ units) with single scat- 

720 tering albedo (Figure 6) implies that shadow hiding is more efficient in high 

721 albedo surfaces. This trend was previously remarked in the opposition slope 
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of asteroids by (iBelskaya and Shevchenkd . 120001 ). It was first interpreted by 
these authors as a decr ease of the absolu t e slop e with albedo, consistent with 



the analytical model of 



Helfenstein et al. 



(119971 ) which predicted that the am- 



725 plitude of the shadow hiding must decrease wit h albedo. However, because the 



726 slope unit in (IBelskaya and Shevchenko 



20001) is magnitude (remark that the 



scale of the mangitude is not reversed for their graph, figure 4, as for the the 
other graphes that show the phase curves, figures 1 and 2), a decreasing slope 
in mag.deg"^ corresponds to an increasing slope in I/F.deg~^ units. As a con- 
sequence, the behavior of the slope as a function of the albedo for the rings, 
satellites and asteroids of the Solar System is consistent and all lead to the 
same idea that shadow hiding is reinforced at high albedo. The correlation be- 
tween slope and albedo seems to be the strongest trend of the opposition effect 
in satellites and rings of the Solar System and the use of more sophisticated 
models is needed to understand them. 



he results from Figure 6 are in agreement with the simulations of ray-tracing. 



737 ( IStankevich et al. 



19991 ). which models shadow hiding in a layer of particles. 
These simulations show that shadow hiding creates a linear part in the phase 
function from 10 to 40 degrees and that the absolute slope of the linear part 
becomes steeper when optical depth increases and the filling factor of the layer 
of particles increases. 

How does albedo relate to the optical depth and the filling factor? Previous 
studies have shown that the albedo and optical depth are highly positively 



corre lated for the rings, (see 



Doyle et al. 



19891 : 



Cooke 



1991 



Pones et al 



19931 ). For satellites, optical depth is effectively infinite; since this removes one 
variable, relating the slope parameter to the nature of the surface is easier for 
satellites than for rings. Thus we must consider two kinds of objects: 
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748 • For rings, where the optical depth is finite, variation of the slope S will be 

749 a subtle effect involving both optical depth and filling factor ; 

750 • For satellites, which have "solid" surfaces, variations in slope are linked to the 

751 filling factor of each surfac e. If the optical d e pth i s invariant for satellites. 



752 according to the model of 



Stankevich et al. 



(119991 ). only variations of the 



753 filling factor can explain differences in the slope S. However, the notion of 

754 filling factor is not well suited for satellites; indeed, a description involving 

755 a topographical roughness is more appropriate. 

756 We noticed that when the optical depth is finite, as for the rings, the ef- 

757 fects of slope are stronger with a high albedo than for high albedo satellites. 

758 Consequently the shadow hiding effect for the ring is more efficient than for 

759 satellites and refiects a difference between the three-dimensional aspect of a 

760 layer of particles in the rings and a planetary regolith. 



761 5 Conclusion 



762 The goal of this paper was to understand the role of the viewing conditions 

763 on the morphological parameters of the opposition surge and the role played 

764 by the single scattering albedo on the morphological parameters. We have use 

765 three methods to fit the data: the first one with a simple morphological model, 

766 secondly by taking more accurately into account the role of the size of the Sun 

767 in the morphological model and thirdly eliminating the role of the size of the 

768 Sin in the morphological model. 

769 The results of this study allow us to highlight several facts related to the 

770 observation and the mechanisms of the opposition effect in the Solar System. 
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'1) The slope of the hnear part is an increasinK function o 



are consistent with those of 



Belskaya and Shevchenkd (120001 ) . for which the 



albedo. Our results 



773 slope of the phase function of asteroids increases when albedo increases. 

774 These results confirm the predictions from simulations of shadow hiding for 

775 the first time. 

775(2) We note that the morphological parameters of the surge (A and HWHM) 

777 are sensitive to the phase angle coverage, specifically to the smallest phase 

778 angles. We have extrapolated observational data points in order to correct 

779 the lack of data near the solar angular radius. However, this method needs 
7B0 to be improved, for example by taking directly into account of the solar 

781 angular radius in the linear-exponential function. We hope that future data 

782 at the smallest phase angles, will confirm the extrapolated data that we use 

783 to perform the "extrapolated linear-exponential" model. 

784(3) The amplitude and the angular width of the opposition surge are linked 

785 to the single scatterin g albedo of the sur faces, as already noted in labo- 



786 ratory measurements (iKaasalaineru . 



20031). Like 



Belskaya and Shevchenko 



787 fj2000l ). we believe that the single scattering albedo is one of the key el- 

788 ements constraining morphological parameters. However before physically 

789 interpreting these results, A and HWHM need to be deconvolved to the "so- 

790 lar size bias" since we have a large dispersion in the relations of A = /(cuq) 

791 and HWHM=/(roo). 

792(4) By deconvolving the phase functions to the Sun's angular size effect, we 

793 showed that A and HWHM are still correlated with the albedo (with better 

794 correlation coefficients). The dependence of A and HWHM are now inde- 

795 pendent of the distance from the Sun, unlike their convolved counterparts. 

796 Indeed, values of A and HWHM from deconvolved phase functions can be 

797 classified into three groups that include a mix of bodies from the inner and 
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798 the outer Solar System. This shows that icy and young surfaces (such as 

799 Europa, lo, Enceladus and Triton) have the smallest amplitudes, whereas 

800 dark and older surfaces (such as the Moon, Phoebe and the C ring) have 

801 the largest amplitudes. 

802(5) It seems that two effects (the "solar size bias" and the "planetary environ- 

803 mental effect"), act together to disperse data taken from different places in 

804 the Solar System. Moreover, with our technique of deconvolution of phase 

805 curves, we see that the "solar size bias" can be removed from A and HWHM, 

806 because unconvolved data have A and HWHM that don't show any trend 

807 with distance from the Sun. These arguments strengthen the conclusions 

808 that the notion of "ecosystem" for a planetary environment can be the key 

809 element determining the opposition effect surge morphology. 

810 Our method cannot directly derive the physical properties obtained from the 

811 models. Firstly, because there is a large set of models and it seemed more 

812 convenient to separate the morphological models from the more physical and 

813 sophisticated ones. Secondly, because the various spectral resolution of our 

814 data set is not appropriate for a majority of physical models which need a fine 

815 spectral resolution (for example, in the coherent backscatter theory, HWHM 

816 is linked to the ratio of the wavelength over the free mean path of photons). In 

817 addition, the coherent backscatter can singificantly polarized the brightness 

818 of a surface (??), so the polarized phase curves can bring crucial and comple- 

819 mentary informations to that of the unpolarized phase curves. Consequently, 

820 more investigations need to be provided for this purpose by using color and 

821 polarized phase curves. 

822 For a future work, which will critically depend on the quality of the observa- 

823 tions, first it would be interesting to study the phase functions of the lead- 
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ing and trailing faces of synchronously rotating satellite in order to test the 
role of the environmental effect more preciselys. Indeed, satellites are sub- 
ject to energy fluxes from electrons, photons and magnetospheric plasma, 
and io n bombardment, wh ich are not the same on the leading and trailing 



828 sides (IBuratti et al. 



1988[ ). Consequently, morphological parameters might 



vary significantly from the leading side to the trailing side for the same satel- 
lites. However, the impo rtant dispersion in trailing side data (see for example 



Kaasalainen et al. 



200 ll ) did not allow us to pursue this comparison. We then 
hope to have in the future that more accurate data of all satellites of the satel- 
lites of the Solar System. 

Secondly, to better understand the role of the "planetary environmental effect," 
a more relevant study would be the comparison of rings with "ringmoons" or 
small satellites which are in the vicinity of the rings. Several examples of such 
a ring/ringmoon system are present in the environment of each giant planet : 



for Jupiter : Metis and Adrastea with the main ring (IShowalter et al.l . 119871 ): 



839 Amalthea and Thebe with the Gossamer ring (IBurns et al. 



1999f) : 



840 • for Saturn : Pan, Dapnh is and Atlas with the outer A ring (ISmith et al 



1981 



Spitale et al. 



20061 ) , Prometheus and Pandora with the F ring ( 



Smith et al. 



198ll ). and Enceladus with the E ring (although Enceladu s is the primary 



843 sourc e of the E ring, it is usually not called a "ringmoon" (jVerbiscer et al. 



20071) ; 



845 • for U ranus : Cordelia and Ophelia with the e ring (IFrench and Nicholson 



19951) ; 



847 • for Neptune : Galatea with the Adams ring (jPorco 



19911 ). 
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848 Unfortunately, opposition phase curves of the small satellites are not actually 

849 available because they require a fine spatial resolution. 
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Fig. 1. Phase curves of a selection of rings and satellites in the Solar System (see 
Table 1 for references). The solid curves correspond to the best fit obtained with 
the linear-by-parts model and the dotted curves to the best linear-exponential fit. 
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Fig. 2. Phase curves of a selection of rings and satellites in the Solar System. The 
solid curves correspond to the best fit obtained with the linear-exponential model 
convolved with the size of the Sun (using extrapolated data below the angular 
size of the Sun, empty symbols) and the dashed curves correspond to the best 
linear-exponential fit (using only original data, filled symbols). The vertical dotted 
lines represent the angular size of the Sun at the observation time (see Table 6 of 
Appendix). 5q 




Fig. 3. Phase curves of a selection of rings and satellites in the Solar System. The 
solid lines correspond to the best fit obtained with the linear-by-parts fit to the 
logarithmic model (in solid curves). 
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Fig. 4. HWHM of the surge for satellites of Jupiter, Saturn, Uranus and Neptune 
derived with : (a) the linear-by-parts model and (b) the linear-exponential model 
convolved with the size of the Sun. 

In (a) dashed line corresponds to a power-law fit to the data which is 
HWHM~ 0.9 X 0.3^^0. 

In (b) dashed line corresponds to a linear fit HWHM~ 0.15 + 0.19ci7o. 
Dotted lines are empirical functions to the boundaries of the data. 



52 



1.5 - 



* 



1.0 1 

0.00 



(b) Linear-by-parts 



(b) Extrapolated Linear-exponential 



< 
O 

i 1.5 h 



0.25 0.50 0.75 

single scattering albedo m 



1.0 

1 .00 0.00 





1 1 1 1 1 1 1 1 

I* 












:* ' 


^'^h 

i: 



• Europa 
Ganymede 

* Callisto 



Phoebe 
■ Portia group 

• Ariel 
•ft Titania 

♦ Oberon 
□ Nereid 
O Triton 



0.25 0.50 0.75 

single scattering albedo m 



1.00 



Fig. 5. Amplitude of the surge derived with the hnear-by-parts model for satellites 

of Jupiter, Saturn, Uranus and Neptune : (a) the linear-by-parts model and (b) the 

linear-exponential model convolved with the size of the Sun. 

In (a) dashed line corresponds to a power fit to the data A~ 1.65 x 0.72^° 

In (b) dashed line corresponds to the fit A~ 1.75 x 0.72™'^. 

Dotted lines are empirical functions to the data boundaries. 
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Fig. 6. Morphological parameter S derived with the linear-by-parts model for : rings 
(a) and satellites (b) of Jupiter, Saturn, Uranus and Neptune. 
In (a) dashed line corresponds to a power-law fit to the data which is 
S~ 0.001 + omw^ 

In (b) dashed line corresponds to a fit to the data which is S~ 0.001 -hO.Olrog. 
Dotted lines are empirical functions to the data boundaries. 
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Fig. 7. Cross comparison between the morphological parameters of the surge derived 
(a) with the linear-exponential model and (b) with the linear-exponential model 
convolved with the size of the Sun (see tables 7 and 6 of Appendix). 
Dashed and dotted ellipses are arbitrary delimitations of the data points (see text). 
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Fig. 9. Variation of the morphological parameters: the amplitude A of the surge (a) 
and the half-width at half-maximum HWHM (b) derived with the linear-by-parts 
model with the distance from Saturn. 
(Distances are taken from 



Murray and Dermott 



2000) 
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Fig. 10. Variation of the morphological parameters: the amplitude A of the surge 
(a) and the half-width at half-maximum HWHM (b) derived with the logarithmic 



model fitted by the linear- 
(Distances are taken from 
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Solid ellipses in (a) and (b) are arbitrary delimitations of the data points (see text). 
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Table 1 

References for opposition phase curves of the selection of Solar System rings and 



satellites. 





Object 


A 
fnm) 


References 




Moon 


~57U 


fWnitaker. 19d9: Koueier. 1\)66) 




Main ring 


* A fin 

~4DU 


( inrooD et ai.. ZUU4) 




iO 


~0 ( u 


I ivicrijwen ei ai.. lyoo) 


piter 


Europa 


~ouu 


fThomuson and Lockwood. 19921 


1-3 


Ganymede 


~DUU 


f Morrison et al.. 19T4; Millis and ThomDson, 1975; Blanco and Catalano, 19741 








\ X llUllilJoUll ciliU. J-jUCiVWUUU.. ±uc/^ } 




u ring 




(l^rencn et ai.. zUU/l 




B ring 


~DOU 


fFranklin and Cook. 19651 




D rmg (^xio i ) 


D ( Z 


(rrencn ex ai.. zuu/ ) 




A ring 


672 


(French et al.. 2007) 


^ 

;3 

4^ 


E ring 


~DOU 


Ir^ane et ai., lyoo. L/arson. iyo4, onowaiter et ai.. iyyil 


03 


Enceladus 




1 veruiscer et ai.. zuuol 




Rlicaj 




I j-/Oiiiiiiti ue ei) ai.. veiuiooei aiiu. veveiK.ci. LaOij \ 








1 1; i dillvnii dilU. V^OUK. A.ai^\ 




Phoebe 


~d5U 


(Bauer et al.. 20061 




IVillgb 




('K'tirVrmnViVn 9001 

I jxai KObciiKa. zuui i 




Portia group 


~500 


fKarkoschka. 2001) 


Uran 


Ariel 


~bUU 


(jDuratti et al.. iyyz: ivarlcoscnica. zUUi) 




i itania 


~dUU 


(iDuratti et al.. lyyz: Karkosciika. zOOi) 




Oberon 


~600 


fBuratti et al.. 1992: Karkoschka. 2001) 




Fraternite 


~500 


fde Pater et al.. 2005: Ferrari and Brahic. 1994) 


Neptun( 


Egalite 
Nereid 


~500 
~570 


fde Pater et al.. 2005: Ferrari and Brahic. 1994) 
fSchaefer and Tourtellotte. 2001) 




Triton 


~400 


('Buratti et al.. 1991) 
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Table 2 

Morphological parameters of opposition phase curves of Solar System objects. The 

unit of HWHM is the degree and the unit of the slope S is I/F.deg~^ 

Linear-by-parts fit Linear-exponential fit 





UDjecx 


A 


TTTT7TTTV /f 
rl W XliVi 


c 


A 

A 


rl W MM 


c 
o 




Moon 


1.27 


1.21 


0.0261 


1.53 


1.98 


0.0017 




Main ring 


1.23 


2.68 


0.0170 


1.34 


2.03 


4.82x10"*^ 


CD 


lo 


1.25 


0.65 


0.0204 


1.14 


0.83 


0.01500 


Jupit 


Europa 


1.11 


0.41 


0.0095 


1.13 


0.31 


0.0054 




Ganymede 


1.31 


0.46 


0.0125 


1.28 


0.69 


0.0040 




Callisto 


1.25 


0.78 


0.0291 


1.50 


2.07 


0.0024 




C ring 


1.61 


0.15 


0.0342 


1.55 


0.09 


0.0030 




B ring 


1.30 


0.23 


0.0312 


1.28 


0.30 


0.0155 




B ring (HST) 


1.38 


0.14 


0.0278 


1.37 


0.09 


0.0238 




A ring 


1.44 


0.13 


0.0297 


1.44 


0.08 


0.0159 


turn 


E ring 


1.22 


0.40 


0.0903 


1.51 


0.86 


1.19x10"'^ 


C» 


EnrelaHns 

J — /XJ.VjV^AC*VA LAO 


1.15 


0.42 


0.0127 


1.20 


0.29 


0.0113 




Rhea 


1.17 


0.41 


0.0142 


1.14 


0.52 


0.0080 




lapetus 


1.22 


0.25 


0.0348 


1.36 


0.22 


0.0030 




Phoebe 


1.17 


0.41 


0.0485 


1.20 


0.38 


0.0003 




Rings 


1.07 


0.21 


0.0274 


1.08 


0.45 


0.0012 




Portia group 


1.80 


0.20 


0.0015 


1.48 


0.13 


0.0015 


Uran 


Ariel 


1.63 


0.38 


0.0178 


1.60 


0.18 


0.0074 




Titania 


1.67 


0.45 


0.0126 


1.82 


0.27 


0.0032 




Oberon 


1.77 


0.37 


0.0174 


1.81 


0.31 


0.0037 


)tune 


Praternite 


2.16 


2.57 


0.0136 


1.86 


1.50 


1.10x10-^ 


Egalite 


1.72 


1.01 


0.0160 


1.77 


0.67 


8.34x10"^ 


CD 


Nereid 


1.56 


0.50 


0.0042 


1.45 


0.34 


0.0099 




Triton 


1.17 


0.10 


0.0124 


1.25 


0.27 


0.0049 
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Table 3 

Direct output parameters of morphological models using opposition phase curves 
and "ideal" opposition phase curves of Solar System objects. The unit of w is the 



degree and the unit of the slope Ig is I/F.deg 





Object 


Logaritlrmic fit 




Linear-expone 


ntial fit 




Exti 


■apolated linear-exponential fit 




ao 


ai 


Ip 


h 


w 


Ib 


Ip 


■ft, 


w 


Is 






Moon 


0.114 


-0.015 


0.048 


0.090 


1.433 


0.0017 


0.049 


0.088 


1.548 


0.0016 




Main ring 


3.9x10"'' 


-5.6x10"" 


1.1x10"^ 


3.2x10"'' 


1.470 


4.8x10"^ 


2.4x10"'= 


3.2x10"'^ 


0.905 


4.6x10"** 


S 


lo 


0.726 


-0.077 


0.100 


0.699 


0.600 


0.0150 


0.313 


0.652 


0.375 


0.0107 


'p. 


Europa 


0.580 


-0.025 


0.076 


0.572 


0.230 


0.0054 


0.086 


0.573 


0.193 


0.0055 




Ganymede 


0.379 


-0.028 


0.100 


0.349 


0.500 


0.0040 


0.114 


0.355 


0.325 


0.0041 




Callisto 


0.182 


-0.025 


0.070 


0.140 


1.500 


0.0024 


0.091 


0.171 


0.248 


0.0051 




C ring 


0.056 


-0.007 


0.033 


0.059 


0.070 


0.0030 


0.031 


0.056 


0.104 


0.0021 




B ring 


0.506 


-0.050 


0.142 


0.502 


0.223 


0.0155 


0.179 


0.511 


0.143 


0.0174 




B ring (HST) 


0.628 


-0.058 


0.249 


0.656 


0.066 


0.0238 


0.213 


0.647 


0.094 


0.0216 




A ring 


0.390 


-0.046 


0.182 


0.409 


0.060 


0.0159 


0.150 


0.401 


0.093 


0.0143 


a 


E ring 


2.1x10"'^ 


-5.2x10"" 


9.3x10"^ 


1.8x10"'^ 


0.623 


1.2x10"'' 


4.8x10"' 


1.2x10"'' 


1.085 


3.8x10"' 


m 


Enceladus 


0.912 


-0.063 


0.180 


0.895 


0.215 


0.0113 


0.251 


0.899 


0.150 


0.0116 




Rhea 


0.602 


-0.055 


0.085 


0.573 


0.379 


0.0080 


0.236 


0.576 


0.139 


0.0082 




lapetus 


0.112 


-0.010 


0.039 


0.110 


0.160 


0.0030 


0.035 


0.116 


0.111 


0.0039 




Plioebe 


0.007 


-0.0001 


0.001 


0.008 


0.276 


0.0003 


0.002 


0.008 


0.067 


0.0004 




Rings 


0.038 


-0.004 


0.003 


0.044 


0.326 


0.0012 


0.014 


0.045 


0.017 


0.0012 


n 


Portia group 


0.058 


-0.006 


0.026 


0.055 


0.100 


0.0015 


0.038 


0.054 


0.273 


0.0010 


a 

d 

S-H 


Ariel 


0.366 


-0.044 


0.213 


0.355 


0.129 


0.0074 


0.217 


0.359 


0.112 


0.0077 


Titania 


0.252 


-0.034 


0.181 


0.221 


0.198 


0.0032 


0.181 


0.221 


0.198 


0.0032 




Oberon 


0.236 


-0.034 


0.168 


0.208 


0.229 


0.0037 


0.181 


0.211 


0.186 


0.0040 


o 


Fraternitc 


0.0005 


-0.0001 


0.00037 


0.00043 


1.082 


1.1x10"^ 


0.00086 


0.00022 


1.1x10"^ 


5.6x10"'^ 


a 


Egalite 


0.0001 


-0.001 


0.00040 


0.00052 


0.489 


8.3x10"'^ 


0.00090 


0.00021 


1.106 


1.5x10"^ 


p. 

o 


Nereid 


0.144 


-0.021 


0.064 


0.140 


0.250 


0.0099 


0.065 


0.175 


0.057 


0.0301 




Triton 


0.584 


-0.035 


0.140 


0.560 


0.200 


0.0049 


0.146 


0.609 


0.049 


0.0075 
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Table 4 

References for the single scattering albedo of Solar System objects. 





Object 




A 
(nm) 


References 




Moon 


0,21 


~500 


fHelfenstein et al.. 1997) 




lo 


0,75 


590 


fMcEwen et al.. 1988) 


>-< 

CD 


Europa 


0,96 


550 


(DominEue and Verbiscer. 19971 


+^ 
■ft 
1-3 


Ganymede 


0,87 


470 


(DominEue and Verbiscer. 1997) 




Callisto 


0,53 


470 


(DominEue and Verbiscer. 1997) 




C ring 


0,16 


672 


fPrench et al.. 2007) 




B ring 


0,83 


672 


fPoulet et al.. 2002) 




B ring (HST) 


0,85 


672 


(French et al.. 2007) 




A ring 


0,79 


672 


(French et al.. 2007) 


Satu 


Enceladus 


0,99 


480 


(Verbiscer and Veverka. 1991) 




Rhea 


0,86 


480 


(Verbiscer and Veverka, 1989) 




lapetus 


0,16 


480 


(Buratti. 1984) 




Phoebe 


0,06 


480 


(Simonelh et al.. 1999) 




Rings 


0,06 


~500 


(Karkoschka. 2001) 




Portia group 


0,09 


~500 


(Karkoschka. 2001) 


Uran 


Ariel 


0,64 


~500 


(Karkoschka. 2001) 


Titania 


0,48 


~475 


(Veverka et al. 1987) 




Oberon 


0,43 


~500 


(Karkoschka. 2001) 




Fraternite 


0,02 


480 


(Ferrari and Brahic, 1994) 


CD 

+^ 
ft 


Egalite 


0,02 


480 


(Ferrari and Brahic, 1994) 


CD 


Nereid 


0,21 


~500 


(Thomas et al.. 1991) 




Triton 


0,97 


500 


(Lee et al., 1992) 
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Table 5 



Morphological parameters of "ideal" opposition phase curves of Solar System ob- 



jects. The unit of HWHM is 


the degree and the unit of the slope 5 


is I/F.deg ^ 




Object 


Linear-by-parts 


and log fits 


Extrapolated linear-exponential fit 




A 


HWHM 




A 


T TT "X 7TT1V /T 

HWHM 


S 




Moon 


2.18 


0.54 


0.0224 


1.55 


2.14 


0.0016 




Main ring 


2.21 


0.55 


0.0230 


1.76 


1.25 


4.55x10-8 


CD 


lo 


1.89 


0.41 


0.0169 


1.47 


0.52 


0.0107 


Jupit 


Europa 


1.34 


0.16 


0.0064 


1.15 


0.26 


0.0055 




Ganymede 


1.62 


0.29 


0.0117 


1.32 


0.45 


0.0041 




Callisto 


2.21 


0.55 


0.0229 


1.53 


0.34 


0.0051 




C ring 


2.17 


0.53 


0.0221 


1.54 


0.14 


0.0021 




B ring 


1.82 


0.38 


0.0155 


1.35 


0.19 


0.0174 




B ring (HST) 


1.77 


0.36 


0.0146 


1.32 


0.13 


0.0216 




A ring 


1.87 


0.40 


0.0165 


1.37 


0.12 


0.0143 


turn 


E ring 


3.38 


0.99 


0.0450 


2.80 


1.50 


3.76x10-^ 


03 


Enceladus 


1 

1 .00 


n 97 


0.0107 


1 97 
l.Z i 


n 9n 


U.Uilu 




Rhea 


1.76 


0.35 


n ni AA 


1.41 


0.19 


0.0082 




lapetus 


1.76 


0.36 


0.0145 


1.30 


0.15 


0.0039 




Phoebe 


1.99 


0.45 


0.0187 


1.33 


0.09 


0.0004 




Rings 


1.98 


0.45 


0.0186 


1.31 


0.02 


0.0012 


m 


Portia group 


1.80 


0.20 


0.0015 


1.71 


0.37 


0.0010 


Uran 


Ariel 


2.13 


0.52 


0.0215 


1.60 


0.15 


0.0077 




Titania 


2.21 


0.55 


0.0230 


1.81 


0.27 


0.0032 




Oberon 


2.41 


0.63 


0.0267 


1.86 


0.25 


0.0040 


)tune 


Fraternite 


2.72 


0.75 


0.0314 


4.99 


1.53 


5.61x10-^ 


Egalite 


2.39 


0.62 


0.0253 


5.34 


1.53 


1.53x10"^ 


CD 


Nereid 


2.27 


0.57 


0.0241 


1.37 


0.08 


0.0301 




Triton 


1.49 


0.23 


0.0092 


1.23 


0.06 


0.0075 
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APPENDIX : 

REFINEMENTS TO THE SUN'S ANGULAR SIZE EFFECT To compute 
the distance from a Solar System object to the Sun at any given date, it is not 
appropriate to use the semi-major axis. We use a series of equations which take into 
consideration the distance between the Sun and the planet at the approximate date. 
The heliocentric radius rp, the distance from the focus of the ellipse (i.e. the Sun) 
to the planet, is given by: 



where E is the eccentric anomaly, a and e are two of the seven orbital elements 
which define an ellipse in space : a is the mean distance, or the value of the semi- 
major axis of the orbit (average Sun to planet distance); e is the eccentricity of the 
ellipse which describes the orbit (dimensionless); i is the inclination (in degrees), 
or angle between the plane of the ecliptic (the plane of the Earth's orbit about the 
Sun) and the plane of the planets orbit; $7 is the longitude of ascending node (in 
degrees), or the position in the orbit where the elliptical path of the planet passes 
through the plane of the ecliptic, from below the plane to above the plane; u) is the 
longitude of perihelion (in degrees), or the position in the orbit where the planet is 
closest to the Sun; A is the mean longitude (in degrees), the position of the planet 
in the orbit; and M is mean anomaly (in degrees). The mean anomaly gives the 
planet's angular position for a circular orbit with radius equal to the semi major 
axis. It is computed directly from the elements using: 



Kepler's second law states that the radius vector of a planet sweeps out equal areas 
in equal times. The planet must speed up and slow down in its orbit. The true 
anomaly / gives the planet's actual angular position in its orbit. It is the angle 
(at the Sun) between perihelion of the orbit and the current location of the planet. 
To obtain its value, first we compute the eccentric anomaly, E, from M and the 
eccentricity e by using the "Kepler equation": 



Tp = a(l — e cos E) 



(10) 



M = A - a) 
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M = E - eshiE 



(12) 



An expansion to order of the solution to the "Kepler equation" is: 



(13) 



(14) 
(15) 
(16) 
(17) 

(18) 
(19) 



where t is the observation time (Table Appendix6) converted in Julian centuries 
since JD 2451545.0, the zero index quantities are the orbital elements at the epoch 



E = M + {e ) sin M + -e^ sin 2M + -e"^ sin 3M 

8 2 8 



Now, to find the orbital elements of a planet at a specific date, we use: 

a = ao + at 
e = eo + et 
i = io + {a /3600)t 
a) = tDo + (i/3600) t 

n = no+ (r2/36oo) t 



A = Ao + A/3600 + 360 Nr t 



of J2000 (JD 2451545.0) and the dot quan tities are the change per ju. 



Murrav and Dermott 



ian century 



2000). When we 



of the orbital elements (values are taken in 
have the heliocentric radius at the given observation time, we can compute the solar 
angular size a© using and tq, the radius of the Sun: 



Oq = arcsm — 



(20) 



The limb darkening function of 



Pierce and Waddelll ( 1961 



) has been used in the past 



to be convolved with a theoretical opposition effect function (iKawata and Irving 



19741 ). Its formula is: 

W{fi') = ax + bx^l' + cx 



1 - /x' • log ( 1 + ^ 



(21) 



where fi' = cos 6' and 6' varies from to the Sun's angular radius Oq. ax, bx and cx 
are coefficients that depend on the wavelength (values are given in table Appendix7) . 
We perform a normalized convolution of the limb darkening function to the linear- 
exponential function P{a) by doing: 



Table 6 

References for the observational parameters needed to compute the angular size of 
the Sun a© . Bold text corresponds to missing data replaced arbitrarily, we assumed 

an hourly time of 18:00:00.0 

Solar Angular Size Observation Time 
Object References 







«0(°) 


UT 






Moon 


0.263 


1968 December 24 


rWhitaker. 1969) 




Main ring 


0.0506 


2000 December 13 


fThrooD et al.. 2004) 




lo 


0.0498 


1977 December 15 


fLockwood et al.. 1980) 


piter 


Europa 


0.0515 


1976 January 1 


fThomDson and Lockwood. 1992) 


1-3 


Ganymede 


0.0499 


1971 may 1 


fBlanco and Catalano. 1974) 




Callisto 


0.0515 


1976 January 1 


fThomoson and Lockwood. 1992) 




C ring 


0.0295 


2005 January 13 


CFrench et al.. 2007) 




B ring 


0.0266 


1959 June 26 


fFranklin and Cook. 1965) 




B ring (HST) 


0.0295 


2005 January 13 


fFrench et al.. 2007) 




A ring 


0.0295 


2005 January 13 


fFrench et al.. 2007) 


,turn 


E ring 


0.0271 


1980 January 1 


fLarson. 1984) 


a: 
CO 


Enceladus 


0.0287 


1997 October 10 


rVerbiscer et al.. 2005) 




Rhea 


0.0268 


1976 January 13 


CLockwood et al.. 1980) 




lapetus 


0.0266 


1972 December 1 


fFrankhn and Cook. 1974) 




Phoebe 


0.0295 


2005 January 13 


f Bauer et al.. 2006) 




Rings 


0.0138 


1997 July 29 


fKarkoschka. 2001) 




Portia group 


0.0138 


1997 July 29 


fKarkoschka. 2001) 


Uran 


Ariel 


0.0138 


1997 July 29 


fKarkoschka. 2001) 




Titania 


0.0138 


1997 July 29 


fKarkoschka. 2001) 




Oberon 


0.0138 


1997 July 29 


fKarkoschka. 2001) 


CD 


Fraternite 


0.00885 


2002 July 27 


fde Pater et al.. 2005) 


4^ 
ft 


Egalite 


0.00885 


2002 July 27 


fde Pater et al.. 2005) 


CD 


Nereid 


0.00889 


1998 June 20 


fSchaefer and Tourtellotte. 2001) 




Triton 


0.00885 


1988 June 20 


fBuratti et al.. 1991) 



Table 7 

Parameters of the limb darkening function of 
wavelengths to the observations. 



Pierce and Waddelll (|l96lh at similar 



Pierce and Waddeh fl9611 







^obs 

(nm) 


\ 

/V 

(nm) 


OA 








Moon 


~570 


560 


0.75079 


0.41593 


-0.54334 




Main ring 


~460 


460 


0.58274 


0.56078 


-0.46772 




Id 


~570 


560 


0.75079 


0.41593 


-0.54334 


piter 


Europa 


~500 


500 


0.68897 


0.47873 


-0.54651 


1-3 


Ganymede 


~600 


600 


0.78074 


0.38427 


-0.53777 




Callisto 


~500 


500 


0.68897 


0.47873 


-0.54651 




C ring 


672 


660 


0.83717 


0.33283 


-0.55400 




B ring 


~650 


640 


0.81999 


0.34918 


-0.55132 




B ring (HST) 


672 


660 


0.83717 


0.33283 


-0.55400 




A ring 


672 


660 


0.83717 


0.33283 


-0.55400 


turn 


E ring 


~650 


640 


0.81999 


0.34918 


-0.55132 


CO 


Enceladus 


439 


440 


0.49375 


0.62584 


-0.38974 




Rhea 


~500 


500 


0.68897 


0.47873 


-0.54651 




lapetus 


~600 


600 


0.78074 


0.38427 


-0.53777 




Phoebe 


~650 


640 


0.81999 


0.34918 


-0.55132 




Rings 


~500 


500 


0.68897 


0.47873 


-0.54651 


in 


Portia group 


~500 


500 


0.68897 


0.47873 


-0.54651 


Uran 


Ariel 


~600 


600 


0.78074 


0.38427 


-0.53777 




Titania 


~600 


600 


0.78074 


0.38427 


-0.53777 




Oberon 


~600 


600 


0.78074 


0.38427 


-0.53777 


)tune 


Fraternite 


~500 


500 


0.68897 


0.47873 


-0.54651 


Egalite 


~500 


500 


0.68897 


0.47873 


-0.54651 


CD 


Nereid 


~570 


560 


0.75079 


0.41593 


-0.54334 




Triton 


~400 


400 


0.16732 


0.84347 


-0.03516 



-6^ 
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■ lo 
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0.00 
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1.00 
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n r 
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